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PROCEDE DE LOCALISATION D UN ELEMENT D'INTERET CONTENU DANS UN OBJET TRID1MENSIONNEL, EN 
PARTICULIER LORS D'UN EX A MEN DE STEREOTAX IE EN MAMMOGRAPHY 



(57) Abstract 

As stereotactic images are digitised, the method consists in the following steps: 
selecting (400) a target pixel in a target region of interest, generating (401). around the 
selected pixel a target window of selected dimensional characteristics and containing the 
said target region of interest, determining (402) in a second image, based on a predetermined 
criterion of selection, a set of pixels, generating (403) around each selected pixel a second 
window of the same dimensional characteristics as the said target window, effecting a 
correlation processing (404) between the grey levels of the pixels of each second window 
and the grey levels of the target window pixels, so as to obtain a correlation value for each 
second window, and identifying the region of interest homologous with the target region 
of interest based on the analysis of the whole set of correlation values obtained, so as to 
minimise the risks of pairing errors between homologous regions of interest. Then the 
element of interest is located from the positions of the two homologous regions. 

(57) Abreg* -a 

Les images stereotax iques etant numensees. on selectionne (400) un pixel ciblc dans 
une region d'intertt cible, on gdnere (401), autour du pixel cible selectionne. une fenetre 
ciblc de caracteristiques dimensionnclles choisics et contenant ladite nfgion d'mterfct cible, 
on determine (402) dans une deuxieme image, selon un critere de selection predetermine, 
un ensemble de pixels, on genere (403) autour de chaque pixel selectionne une deuxieme 
fenetrc de memes caracteristiques dimensionnelles que ladite fenetre cible, on effectue un 
traitcment de correlation (404) entre les niveaux de gris des pixels de chaque deuxieme 
fenetre et les niveaux de gris des pixels de la fenetre cible. de facon a obtenir une valeur de 
correlation pour chaque deuxieme fenetre, et on identifie la region d'mtertt homologue de 
la region d'interet cible a partir de l'analyse de l'ensemble des valeurs de correlation ainsi 
obtcnues. de facon a minimiscr les risques d'erreur d'appariement entre les regions d intettt 
homologues. On localise ensutte reiement d'interet a partir des positions des deux regions homologies. 
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Procede de localisation d'un element d'interet contenu dans un 
objet tridimensionnel, en particulier lors d'un examen de 
stereotaxic en mammographie. 

L'invention concerne la localisation d'un Element d'interet contenu 
dans un objet tridimensionnel a partir des positions de regions 
d'interet homologues correspondant audit element d'interet et 
apparaissant dans un ensemble d'images stereotaxiques dudit objet. 

5 Bien que l'invention s'applique a tout ensemble d'images 

stereotaxiques d'un objet tridimensionnel quelconque, elle est 
particulierement utile dans le domaine medical et notamment lors d'un 
examen de stereotaxic en mammographie pour la locabsation de 
microcalcifications dans un sein. 

10 Un examen de st£r£otaxie en g6n£ral permet d'acc£der avec une 

grande precision, de l'ordre du millimetre, a un point donn6 d'un objet 
tridimensionnel, a partir de deux projections bidimensionnelles de cet 
objet acquises par exemple selon deux incidences angulaires opposees 
de part et d'autre de la normale au plan du recepteur d'images d^livrant 

15 lesdites projections. 

A partir de la connaissance, d'une part, de la position de la 
projection de ce point dans le plan de chaque image stereotaxique 
obtenue, et, d'autre part, de la geomdtrie de l'appareil de prises de 
vues stereotaxiques ayant conduit a la construction des deux images, il 

20 est possible de calculer la position exacte spatiale de ce point dans 

1'objet tridimensionnel, par calcul trigonomdtrique. 
En resume, un examen stfireotaxique n6cessite : 

de faire au moins deux images de l'objet k diff£rentes angulations, 
de connaitre parfaitement la g€om6trie d^cquisition du systfcme 
25 st6reotaxique, et 

de pouvoir localiser sur les differentes images stereotaxiques 
obtenues, le lieu de la projection de l'element d'interet choisi. 
Les deux premiers points ne posent pas de probleme majeur. 
Actuellement, en mammographie, l'examen de st6r£otaxie est 
30 effectue & l'aide d'un mammographe equip6 d un dispositif de prises de 

vues stereotaxiques. Le mammographe comprend un tube h. rayons X, 
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situe a l'extremite d'un premier bras qui est mobile autour d un axe el 
emettant un rayonnement X vers un recepteur situe a l'extremite d'un 
autre bras. Entre le tube et le recepteur sont disposes, d'une part, un 
plateau porte-sein, ou support patient, et, d'autre part, un plateau de 
compression assurant le maintien en place du sein lors des 
mammographies. Le recepteur d'images peut etre un recepteur 
numerique, tel qu'une camera CCD par exemple, delivrant des images 
ster6otaxiques numerisees. Les prises de vues st6r6otaxiques 
necessitent la rotation du tube a rayons X autour des plateaux porte- 
sein et de compression, selon deux orientations successives opposees 
de part et d'autre de sa position initiale perpendiculaire au plan du 
r6cepteur d'images. 

A partir des images stereotaxiques numerisees obtenues, la 
localisation de l'element d'interet sur les differentes images est une 
operation delicate, en particulier en mammographie, puisqu'elle fait 
appel a une mise en correspondence, ou appariement ("matching" en 
langue anglaise) de regions d'interet homologues correspondant, dans 
deux images ster6otaxiques differentes, a un meme element d'interet 
situ€ dans le sein. 

Ainsi, en mammographie, on cherche a apparier sur les differentes 
images stereotaxiques obtenues, les regions d'interet qui correspondent 
a une meme microcalcification. 

Or, lappariement des microcalcifications projetees est un 
probleme difficile car les microcalcifications ne se ressemblent pas 
d'une image a l'autre. Leur forme et leur contraste peuvent changer 
ainsi que leur agencement dans l'espace. Elles peuvent se superposer 
avec des zones fibreuses ou avec d'autres microcalcifications. 

Jusqu'a present, 1'appariement etait effectu6 entierement 
manuellement par le radiologue. Ce dernier choisissait quelles etaient 
les regions d'intergt en correspondance sur les deux images 
stereotaxiques angulees, puis rentrait les coordonnees 
bidimensionnelles de ces regions dans un calculateur qui donnait la 
position spatiale de la microcalcification, c'est-a-dire les coordonnees 
tridimensionnelles de cette derniere. A partir de ces coordonndes, le 
radiologue pouvait ajuster par exemple la position d'un porte-aiguille 
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de fagon a ce que 1'aiguille vienne h l'aplomb de la microcalcification 
pour perrnettre sa ponction par exemple. 

Or, cet appariemment entierement manuel s'avere etre une 6tape 
longue et parfois imprecise. 

L'invention vise a apporter une solution & ce probleme. 

L'invention a pour but de r£duire sensiblement la duree d'un 
examen st£reotaxique et d'obtenir une meilleure precision 
d'appariement. 

D'une fa9on tres generale, selon l'invention, on d6finit une mesure 
de similarity permettant d'evaluer la ressemblance entre une region 
d'interet (microcalcification par exemple) s£lectionnee (manuellement 
par le radiologue par exemple, ou bien automatiquement par un 
procede de selection automatique) sur une des images ster£otaxiques 
angulees et une region d'interet d'une deuxieme image st£r£otaxique. 

A partir de cette mesure de similarity, on selectionne, pour obtenir 
la region d'interet homologue, un certain nombre de "candidats" 
appartenant a une zone de recherche et ayant une forte ressemblance 
avec la microcalcification de depart. 

Ensuite, en utilisant les principes de la st£r6otaxie, on calcule les 
coordonnees tridimensionnelles de ces candidats, puis on les reprojette 
sur l'image de centrage. On calcule a nouveau la similarity entre ces 
points reprojetes et la microcalcification initialement choisie. En 
effet, si le candidat correspond effectivement a la microcalcification 
de depart, alors sa projection calculee correspondra effectivement a 
l'image de cette microcalcification sur 1'image de centrage. 

Enfin, a l'aide de ces deux valeurs de similarite calculees pour 
chaque candidat, on propose le meilleur candidat. 

En thdorie, on pourrait s'abstenir de 1'etape de reprojection 
utilisant l'image de centrage et se contenter de l'appariement tel qu*il 
vient d'etre expose pour selectionner la region d'interet homologue de 
celle selectionnee sur l'image cible. Cependant, il est nettement 
preferable d'utiliser la reprojection des candidats sur l'image de 
centrage pour verifier la premiere selection effectuee, notamment dans 
le domaine de la mammographie ou les microcalcifications peuvent 
changer de forme et de contraste d'une image h l'autre. 
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En d'autres termes, l'invention propose tout d'abord un precede de 
localisation d'un element d'interet contenu dans un objet 
tridimensionnel a partir des positions de regions d'interets homologues 
correspondant audit element d'interet et apparaissant dans un ensemble 
d'images stereotaxiques dudit objet. Ce procede comporte une etape de 
selection dans une premiere image stenSotaxique, ou image cible, d'une 
premiere region d'interet, en particulier une microcalcification, dite 
region d'interet cible. Le proc6d6 comporte egalement un appariement 
de ladite premiere region avec une deuxieme region d'interet 
homologue de la premiere et apparaissant dans une deuxieme image 
st6r6otaxique. 

Selon une caracteristique generale de l'invention, les images 
stereotaxiques etant num6ris6es, on seiectionne dans 1'etape de 
selection un pixel cible dans ladite region d'interet cible. Dans l'etape 
dappariement, on genere autour du pixel cible seiectionne, une fenetre 
cible de caracteristiques dimensionnelles choisies et contenant ladite 
region d'interet cible. On determine alors dans la deuxieme image, 
selon un critere de selection predetermine, un ensemble de pixels et on 
genere autour de chaque pixel selectionne une deuxieme fenetre de 
memes caracteristiques dimensionnelles que ladite fenetre cible. On 
effectue un traitement de correlation entre les niveaux de gris des 
pixels de chaque deuxieme fenetre et les niveaux de gris des pixels de 
la fenetre cible de fagon a obtenir une valeur de correlation pour 
chaque deuxieme fenetre. On identifie alors la region d'interet 
homologue de la region d'interet cible a partir de l'analyse de 
Tensemble des valeurs de correlation ainsi obtenues. On minimise 
ainsi les risques d'erreur d'appariement entre des regions d'interet 
homologues. 

Selon un mode de mise en oeuvre, l'analyse des valeurs de 
correlation obtenues comporte la selection d'un nombre choisi de 
maxima ou de minima de correlation, la selection de la region 
homologue s'effectuant parmi celles dont la valeur de correlation 
associee est I'un de ces maxima ou minima de correlation. 

En fait suivant que l'image est n normale M ou inversee (une image 
"normale" etant definie comme presentant des niveaux de gris et des 
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contrastes qui sont semblables a ceux habituellement visibles sur des 
cliches de film) on utilisera respectivement les maxima ou les minima. 

On peut egalement prevoir une determination de la dynamique des 
maxima ou minima sdlectionnds, et la comparaison des valeurs de 
dynamique obtenues & un seuil. 

Le traitement de correlation utilise peut comporter un traitement 
de correlation normee, centr6e ou non, ou bien un traitement de 
difference norm6e. 

En variante a Tappariement evoque ci-dessus, l'invention propose 
egalement un appariement base sur 1'eiaboration de caracteristiques 
pr6determin6es de toute region d'inter&t, en particulier des 
caracteristiques physiques comme Tdlongation, le contraste, le 
gradient, puis une comparaison entre ces listes de caracteristiques 
associees a differentes regions d'interSts. 

Plus precis6ment, selon une caracteristique generale de 
l'invention, on selectionne dans l'etape de selection, un pixel cible 
dans ladite region d'interet cible et, dans l'6tape d'appariement, on 
genere, autour du pixel cible selectionne, une fenetre cible de 
caracteristiques dimensionnelles choisies et contenant ladite region 
d'interet cible. On determine pour ladite region d'inter£t cible, a partir 
des valeurs de niveau de gris des pixels de ladite fenfetre cible, un 
premier ensemble de valeurs numeriques respectivement 
representatives de caracteristiques predeterminees caracterisant toute 
region d'interet (par exemple des caracteristiques de forme de 
contraste, de gradient ... ) et on stocke ce premier ensemble de valeurs 
numeriques. On determine dans la deuxieme image, selon un critere de 
selection predetermine, un ensemble de pixels appartenant 
respectivement a des regions d^nterets dites candidates apparaissant 
dans la deuxieme image et on genere autour de chacun de ces pixels 
seiectionnes une deuxieme fenfetre contenant ladite region d f interet 
candidate correspondante. On determine, pour chaque region d'interet 
candidate a partir des valeurs de niveau de gris des pixels de la 
deuxieme fenetre associee, un deuxieme ensemble de valeurs 
numeriques respectivement representatives desdites caracteristiques 
predeterminees et on stocke ce deuxieme ensemble de valeurs 
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numeriques. On identifie alors la region d'interet homologue de la 
region d'interet cible a partir d'un traitement de comparaison entre le 
premier ensemble de valeurs numeriques et chaque deuxieme ensemble 
de valeurs numeriques. On minimise egalement ainsi les risques 
d erreur d'appariement entre des regions homologues. 

Le traitement de comparaison peut comporter un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

Dans les deux variantes d'appariement qui viennent d'etre 
evoqu6es, et afin de minimiser encore le temps de calcul, on determine 
dans la deuxieme image stereotaxique, une zone epipolaire contenant 
au moins le segment epipolaire relatif au pixel cible. Les pixels de la 
zone epipolaire contiennent alors les pixels selectionnes. Pour reduire 
encore le temps de calcul, au lieu d'effectuer l'appariement sur tous les 
pixels de la zone epipolaire determin£e, on seiectionne 
avantageusement les pixels de la deuxieme image parmi ceux 
presentant des maxima ou des minima de niveau de gris dont la 
dynamique est superieure & un seuil predetermine. 

En outre, pour s'affranchir des probl&mes de fond d'image, et pour 
minimiser encore les risques derreurs d'appariement entre deux 
regions d'interets homologues, il est particulierement avantageux que 
l'appariement selon l'invention comporte prealablement un filtrage des 
images stereotaxiques, par exemple du type dit "chapeau haut de 
forme". 

L'utilisation de 1'etape de reprojection des elements d'interdts 
candidats dans Timage de centrage est independante du type 
d'appariement utilise prealablement. 

Ainsi, selon une caracteristique generale de l'invention, le procede 
de localisation du type de celui evoque ci-avant comprend une 
selection dans une premiere image stereotaxique numerisee d'une 
premiere region d'interet dite cible, une selection, dans une deuxieme 
image stereotaxique numerisee, a partir d'un premier appariement 
automatique, d'au moins une deuxieme region d'interet susceptible 
d'etre homologue de la region cible, la determination de la position 
spatiale d'un pixel candidat d'un element d'interet candidat 
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correspondant a ces deux regions d'interets, la determination de la 
position dans une troisieme image st£reotaxique, d'un pixel projete 
correspondant a la projection dans cette troisieme image stereotaxique 
du pixel candidat, et un deuxieme appariement automatique entre la 

5 region int6r6t cible et un voisinage du pixel projete definissant une 

region d'int6ret projet^. 

En pratique, on selectionne avantageusement dans la deuxieme 
image ster6otaxique, h partir dudit premier appariement automatique, 
un ensemble de deuxiemes regions susceptibles d'etre homologues de 

10 la region cible et on determine la position spatiale de chaque pixel 

candidat correspondant a chaque paire de regions d'interets form6e par 
la region cible et Tune des deuxiemes regions. On determine 
dgalement la position dans la troisieme image de chaque pixel projete 
correspondant et on effectue le deuxieme appariement automatique 

15 entre la region cible et un voisinage de chaque pixel projete. 

On peut ordonner Tensemble des deuxiemes regions obtenues dans 
J'ordre decroissant de leur probability a etre la region homologue de la 
region cible. On peut egalement ordonner 1'ensemble des regions 
projetees obtenues dans Tordre decroissant de leur probability a 6tre la 

20 region homologue de la region cible. On selectionne alors 

avantageusement la region homologue de la region cible comme etant 
celle dont le produit des rangs respectifs dans les deux 
ordonnancements est mininum. 

Le premier appariement automatique peut etre l'un des deux 

25 appariements qui ont ete evoqu£s ci-avant, c'est-a-dire a base d'un 

traitement de correlation ou d'un traitement de comparaison. 

De meme, le deuxieme appariement automatique peut cornporter 
un traitement de correlation tel que defini ci-avant, entre le voisinage 
de chaque pixel projetd et une fenStre cible contenant la region 

30 d'interet cible. Ce deuxieme appariement automatique peut etre 

egalement un traitement de comparaison tel que d£fini ci-avant, entre 
la region d'interet cible et chaque region d'int£ret projetee. 

D'un point de vue materiel, 1'invention s'analyse comme un 
dispositif de localisation d'un element d'interet contenu dans un objet 

35 tridimensionnel a partir des positions de regions d'interets homologues 



WO 98/01826 



8 



PCT/FR97/01232 



correspondant audit Element d'intdret et apparaissarit dans un ensemble 
d'images stereotaxiques dudit objet. Ce dispositif comporte un 
appareillage de prises de vues stereotaxiques 6quipe d'un support pour 
ledit objet et d'un recepteur damage tel qu'une camera CCD capable de 
fournir des images stereotaxiques numeris6es. II est egalement prevu 
des moyens de visualisation de ces images, par exemple un ecran 
video, tel que celui d'un microordinateur. Des moyens de selection, 
par exemple une "souris" d'un microordinateur permettant de deplacer 
un marqueur sur I'^cran et de "cliquer" ce marqueur sur un pixel cible 
d'une region d'interet, permettent la selection du pixel cible. II est 
enfin prevu un microprocesseur et des mdmoires associees" pour 
incorporer sous forme logicielle des moyens permettant de r^aliser 
fonctionnellement toutes les Stapes des proced6s qui ont 6x6 decrites 
ci-avant. 

D'autres avantages et caracteristiques de l'invention apparaitront a 
1'examen de modes de mise en oeuvre nullement limitatifs, et des 
dessins annexes, sur lesquels : 

la figure 1 est une vue schematique de trois images stereotaxiques 
d'un objet tridimensionnel a partir d'un dispositif de prises de vues 
stereotaxiques; 

la figure 2 illustre tres schematiquement la notion de segment 
6pipolaire; 

la figure 3 illustre une zone epipolaire choisie conform6ment a 
l'invention pour effectuer un appariement; 

la figure 4 est un organigramrne general d'un procede selon 
l'invention; 

la figure 5 illustre plus particulifcrement une etape d'une variante 
dii proced6 utilisant un traitement de correlation; et 

la figure 6 illustre schematiquement plus particulierement un 
procede selon l'invention utilisant une reprojection sur Tirnage de 
centrage. 

Un examen stereotaxique, en particulier en mammographie, se 
compose d'une serie de trois expositions d un objet tridimensionnel 1 
(figure 1), par exemple un sein, reposant sur un support 3 et comprime 
par un plateau de compression 2, a l'aide d'un tube a rayons X 
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occupant respectivement trois positions differentes 6, 7 et 8. En 
pratique, on procede a une exposition selon un angle de 0° et a deux 
expositions angulees selon deux angles opposes egaux, en pratique 
±15°. 

Un recepteur CCD 4 permet d'obtenir ainsi trois images 
st6reotaxiques numerisees 10, 11 et 12. Selon une convention 
habituellement utilisee, I'image 10 est I'image de droite, tandis que 
I'image 12 est I'image de gauche et que 1'image 11 est 1'image de 
centrage. 

Un element d'int^ret 5, par exemple une microcalcification, 
contenu dans l'objet tridimensionnel 1, fournit sur chacune des images 
10, 11 et 12 des regions d'interets respectivement r6f6rencees 5d, 5c et 
5g. L'image de centrage sert en particulier a verifier que la 
microcalcification sera accessible durant l'acte medical eventuel 
ulterieur. 

Les coordonndes tridimensionnelles de la microcalcification 5 
peuvent etre obtenues, de fa?on classique et bien connue de l'homme 
du metier, par un calcul simple de trigonomdtrie a partir de la 
connaissance des coordonnees bidimensionnelles des deux projections 
5d et 5g, par exemple. Ces coordonnees sont exprimees par rapport a 
une origine materialisee par le centre dune croix 9 gravee sur le 
plateau supportant le sein. Et on repere sur chaque image le centre de 
cette croix de reference pour initialiser geom6triquement le systeme. 

L'homme du mdtier sait que, par definition, le segment dpipolaire 
50 (figure 2) d'un point 5d d'une image 9 d'angulation 91, sur une 
image 12 d'angulation 62, correspond a l'ensemble des points de la 
seconde image 12 pouvant avoir le point 5d comme ant6c6dent. 
Autrement dit, c'est la projection conique par la source d'angulation 
G2, du segment de droite 52 joignant ledit point a la source 
d'angulation 01. 

Dans le cas ou les images stereotaxiques ont 6t6 obtenues avec des 
angulations telles que les positions du tube a rayons X se situent dans 
un plan parallfele aux lignes des images, ce segment dpipolaire 50 est 
un segment de droite qui se trouve etre parallele aux lignes de I'image. 
Cependant, si Ton n'a pas la connaissance de lepaisseur de l'objet 
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tridimensionnel, on ne peut pas determiner la position des extn£mit6s 
de ce segment epipolaire. Ce segment epipolaire devient alors une 
droite epipolaire qui, dans la geom£trie d'acquisition particuliere 
evoquee ici, se situe sur la meme ligne d'image 12 que le point 
d'interet 5d sur l'image 9. 

En outre, lorsqu'on utilise un detecteur qui pr6sente des 
distorsions g6om6triques, la droite epipolaire se transforme en une 
courbe Epipolaire 53 (figure 3). On peut determiner les coordonnees de 
l'ensemble des points de la courbe epipolaire sur l'image 12 par une 
calibration geometrique de la camera CCD. II serait alors possible de 
n'effectuer l'appariement selon 1'invention que sur l'ensemble de ces 
points. 

Cependant, si pour des raisons de prEcision de la calibration 
gEomEtrique, on ne peut determiner qu'approximativement les 
coordonnees des points de la courbe Epipolaire, on choisira alors 
d'effectuer l'appariement selon l'invention sur un ruban epipolaire 51 
qui encadre la courbe Epipolaire. En pratique, on dEfinira une zone 
epipolaire, ou zone de recherche, dElimitEe par deux lignes parallEles 
passant respectivement par deux points augments de part et d'autre du 
point le plus haut et le plus bas de l'epipolaire 53 d'une marge 
predeterminee. Generalement, on opte pour une marge de 3 a 5 pixels. 

On va maintenant decrire plus en detail, en se rEfErant plus 
particulierement aux figures 4 et 5, deux types d'appariement selon 
Tinvention. 

La premiere Etape consiste en la selection (Etape 400) d'un pixel 
cible 60d au sein d'une region d'interet cible 5d presente par exemple 
dans l'image de droite 10, dite image cible. En pratique, cette 
selection peut s'effectuer manuellement en cliquant un point de la 
region d'interet de l'image concernEe apparaissant sur un Ecran d'un 
microordinateur. Ceci etant, on peut prevoir egalement un procEdE de 
selection automatique de pixel. Ce pixel 60d peut Etre centre ou non 
dans la region cible. 

On genere ensuite (Etape 401) une fenetre cible 61 d contenant 
ladite region d'intEret cible. Cette fenEtre cible peut etre gEnEree 
automatiquement apres avoir appliquE sur le pixel cible 60d un moyen 
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logiciel classique de croissance de region. En variante, on peut 
egalement prSvoir, selon le domaine d'application de Tinvention, de 
choisir des fenetres de taille predefinies. Dans tous les cas, on 
obtiendra generalement une fenetre cible rectangulaire ou carree dont 

5 les cotes presentent des dimensions de 1'ordre d'une dizaine a quelques 

dizaines de pixels. 

Les niveaux de gris de chacun des pixels de cette fenetre cible 
sont ensuite stockes dans une memoire associee a un microprocesseur 
qui incorpore de fafon logicielle tous les moyens fonctionnels de mise 

10 en oeuvre du procddd decrit. 

A partir de 1&, on va s61ectionner dans une etape 402 un ensemble 
de pixels dans une deuxieme image ster6otaxique, par exemple l'image 
de gauche. Cet ensemble de pixels selectionnds va constituer une zone 
de recherche 51 a Tinterieur de laquelle on va s61ectionner une ou 

15 plusieurs regions d'interets candidates susceptibles d'etre homologues 

de la region d'interet cible. 

Comme on l'a explique ci-avant et dans le but de reduire la duree 
d'appariement, cette zone de recherche 51 sera par exemple la zone 
epipolaire 51 ddfinie ci-avant. Plus precisement, autour de chacun de 

20 ces pixels selectionnes dans la zone de recherche, on g£nere (6tape 

403) une deuxieme fenetre 61 g de meme caract€ristique 
dimensionnelle que la fenetre cible. On stocke les niveaux de gris des 
pixels de chaque deuxieme fenetre et on procede (etape 404) a un 
traitement de correlation entre les niveaux de gris des pixels de chaque 

25 deuxieme fenetre 61 g et les niveaux de gris des pixels de la fenStre 

cible 6 Id. 

Le traitement de correlation permet done d'effectuer une mesure 
de similarite entre la fenetre cible et chacune des deuxiemes fenetres 
de la deuxieme image de fa?on & obtenir une valeur de similarity 
30 (valeur de correlation entre ces deux fenfetres). 

C'est a partir de l'analyse de ces valeurs de correlation que Ton va 
proceder a I'appariement de la r6gion d'interet cible avec sa region 
d'interet homologue. 

En ce qui conceme le traitement de correlation (etape 404), on 
peut uliliser une correlation norm6e. Plus precisement, chaque valeur 
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de correlation d'une deuxieme fenetre sera alors donn6e par la 
formule : 

N 

I (k)J(k) _^ _^ 

*=i I . J 

\/5 i 2 c<)2 J 2 « iml iun 

ou N est le nombre de pixels d'une fenetre, k l'indice courant d'un 
pixel, I(k) le niveau de gris du pixel k de la fenetre cible et J(k) le 
niveau de gris du pixel k de la deuxieme fenetre. 

On peut egalement utiliser une difference normee permettant 
d'obtenir la valeur de correlation par la forrriule : 

n"]T- J^ll 2 

En th6orie, la region d'interet homologue sera celle contenue dans 
la deuxieme fenetre qui pnSsente la valeur de correlation maximale. 
Ceci etant, en pratique, plusieurs maxima de correlation peuvent etre 
trouvds avec une intensity elevee. On selectionne alors (etape 405), 
parmi I'ensemble des valeurs de correlation, un ensemble restreint de 
maxima de correlation selon un critere predetermine comme l'intensite 
du maximum ou comme le nombre de maxima dans la liste des valeurs 
de correlation rangees par ordre d'intensite decroissante. 

On peut egalement utiliser comme critere suppiementaire de 
selection une determination de la dynamique des maxima seiectionnes 
et la comparaison des valeurs de dynamique obtenues a un seuil 
predetermine. 

La determination de la dynamique de maxima est parfaitement 
connue de l'homme du metier. Cependant, celui-ci pourra 
eventuellement se reporter a la demande de brevet fran9ais n° 
91 15308 pour plus de details concernant une telle determination de 
dynamique. 

II convient de noter ici que le procede selon 1'invention permet de 
determiner et d'afficher sur recran du microordinateur, par exemple 
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par un clignotement, la region d'interet consid6r6e comme 6tant la 
region homologue de la region cible, et ce avec le minimum de risque 
d'erreur. 

Ceci etant, en pratique, le precede selon 1'invention pourra 

5 selectionner un nombre restreint de regions dinterets candidates 

susceptibles d'etre hornologues de la region cible. Le radiologue aura 
alors la possibility de selectionner manuellement celle qu'il considfcre 
comme etant effectivement la region homologue. 

Ainsi qu'on Ta 6voque ci-avant, on peut utiliser dans le traitement 

10 d'appariement, un traitement de comparaison specifique au lieu d'un 

traitement de correlation. 

Ainsi, comme illustnS plus prdcisement sur la figure 4, on 
determine (etape 409) pour la region cible contenue a l'interieur de la 
fen€tre cible 6 Id, et a partir des niveaux de gris des pixels contenus 

15 dans cette fenetre cible, une liste de caracteristiques permettant de 

caractenser en general une region dint€r6t. Ces caracteristiques 
peuvent comporter par exemple une valeur representative de 
1'eiongation de la region d'interet, une valeur representative de sa 
largeur moyenne, une largeur representative de son niveau de gris 

20 moyen, ou de son gradient moyen d'intensite. 

On effectue la meme operation (etape 410) pour les regions 
dinterets candidates contenues dans des deuxiemes fenetres 61g 
definies a partir de pixels selectionnes dans la zone de recherche. A 
cet egard, afin de limiter le temps de traitement, les pixels 

25 selectionnes pour la determination des deuxiemes fenetres englobant 

les regions candidates, sont choisis parmi les pixels de la zone de 
recherche presentant des maxima de niveaux de gris dont la dynamique 
est superieure a un seuil predetermine. 

A partir de ces differentes listes de caracteristiques ainsi 

30 obtenues, on effectue un traitement de comparaison 411 entre la liste 

de caracteristiques affectees a la region cible et chacune des listes de 
caracteristiques affectees & chacune des regions candidates. 

En pratique, bien qu'on puisse envisager un autre type de 
traitement de comparaison, on utilisera de preference un traitement de 

35 comparaison prevoyant la minimisation d'une distance, par exemple la 
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distance eudidienne, entre la liste des caracteristiques de la r6gion 
cible, consid6r6e alors comme un vecteur, et une liste de 
caracteristiques d'une region candidate, consid6ree egalement comme 
un vecteur. 

On pourrait Egalement prevoir, a titre d'exemple, de determiner la 
valeur moyenne des valeurs absolues des differences entre les valeurs 
homologues de deux listes de caracteristiques, et de rechercher celle 
des valeurs moyennes qui est la plus faible pour en deduire la region 
d'interet homologue. 

La selection de la region homologue pr^vue dans I'etape 405 
s'obtient done par cette variante en minimisant une distance entre des 
listes de caracteristiques plutot qu'en maximisant une correlation 
normee ou une difference norm6e. 

L'invention prevoit encore certaines ameliorations. 

Ainsi, dans la variante de l'appariement prevoyant un traitement 
de correlation, on peut limiter la zone de recherche aux pixels de la 
zone epipolaire presentant des maxima de niveaux de gris dont la 
dynamique est superieure & un seuil predetermine. 

En outre, il est preferable, notamment pour s'affranchir des 
problemes de fond d'image, d'effectuer un filtrage prealable des 
images stereotaxiques du type dit "chapeau haut de forme" (ou 
transformation top-hat en langue anglaise). Ce genre de filtrage est 
connu de I'homme du metier et celui-ci pourra se referer pour plus de 
details a Tarticle de J. SERRA Image Analysis and Mathematical 
Morphology, Vol. 2, Academic Press 1988. 

Ce genre de filtrage permet d'extraire d'une image numerique les 
zones claires et etroites et permet de s'affranchir du fond. 

Les etapes du precede selon l'invention qui vont maintenant etre 
decrites, et prevoyant notamment une reprojection sur 1'image de 
centrage, s'apparentent ainsi notamment comme etant une verification 
des resultats de l'appariement precedent entre l'image cible et la 
deuxieme image, et permettent ainsi de minimiser encore les risques 
derreurs d'appariement. 

Ceci etant, bien que la suite du procede qui va maintenant 6tre 
decrite en se referant plus particulierement aux figures 4 et 6, pr6voit 
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la possibility d'utiliser des traitements de correlation ou de 
comparaison du type de ceux qui viennent d'etre decrits, ce proc£de 
selon l'invention utilisant une reprojection sur l'image de centrage, est 
independant de l'appariement prealable utilise entre l'image cible et la 

5 deuxieme image. 

On suppose maintenant (figure 6) que le procede selon l'invention 
a permis de selectionner sur l'image de gauche 12 deux regions 
d interets candidates 5g et 20g susceptibles d'etre la region homologue 
de la region cible 5d apparaissant sur l'image cible 10. 

10 A partir du pixel cible de la region d'interet cible et du pixel 

candidat ayant donn6 lieu a la selection de la region d'interet 
candidate 5g, on determine (etape 406), par un calcul trigonometrique 
de st6reotaxie, et connaissant la geometrie du dispositif de stereotaxic, 
les coordonnees tridimensionnelles du point correspondant de l'objet 

15 d'interet candidat correspondant a ces deux regions 5d et 5g. Puis, par 

un calcul trigonometrique de stereotaxic inverse, on determine les 
coordonnees dans l'image de centrage du pixel projete correspondant a 
la projection du pixel de l'dlement candidat dans cette image de 
centrage. 

20 On remarque done que, dans le cas present, puisque la region 

candidate 5g est effectivement la region homologue de la region 5d, la 
position du pixel projete dans l'image de centrage 11 correspond 
effectivement a une region d'interet projete 5p qui est effectivement la 
projection dans cette image de centrage de reiement d'int£r£t 5. 

25 Par contrc, par ces memes calculs, la region d'interet candidate 

20g va conduire a l'obtention d'une position projetee 20p dans l'image 
de centrage qui, dans le cas present, ne se superposera h aucune region 
d'interet physiquement presente dans cette image de centrage. 

Aussi, un traitement de correlation entre la fenetre cible contenant 

30 la region cible 5d et un voisinage du pixel projete contenant la region 

d'interet projetee 5p va donner une valeur de correlation superieure a 
la valeur de correlation obtenue entre la fenetre cible 5d et un 
voisinage du pixel projete 20p. II en aurait ete de meme, meme si le 
voisinage du point projete 20p avait correspondu ii une micro- 

35 calcification projetee autre que celle correspondant effectivement a la 
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microcalcification 5. 

Ce proc£d£ permet done de selectionner parmi les deux candidats 
5g et 20g la region 5g qui est effectivement la region homologue de la 
region 5d. 

Bien entendu, la conclusion aurait 6i6 la meme en utilisant un 
traitement de comparison 412 au lieu d'un traitement de correlation 
407. 

Utilisation d'un voisinage du pixel projet6, typiquement 
quelques pixels autour du point projetd, peiroet de prendre en compte 
les erreurs eventuelles dues a un bouge de l'objet entre Imposition a 
0°et les expositions h ±15°. 

D'une fa9on plus generate, on peut ordonner l'ensemble des 
regions candidates obtenues dans la deuxi&me image dans l'ordre 
decroissant de leur probability a etre la region homologue de la region 
cible. Ainsi, la region ayant par exemple donn6 le maximum de 
conflation avec la dynamique la plus forte sera affectee du rang 1 et 
ainsi de suite par ordre decroissant. On peut effectuer ce meme 
ordonnancement pour l'ensemble des regions projetees obtenues sur 
Timage de centrage. On s61ectionnera alors la region homologue de la 
region cible comme etant celle dont le produit des rangs respectifs 
dans les deux ordonnancements est minimum. 
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REVENDICATIONS 

1. Proc6d6 de localisation d'un element d'interet contenu dans un 
objet tridimensional a partir des positions de regions d'interfets 
homologues correspondant audit element d*int6rSt et apparaissant dans 
un ensemble d'images stereotaxiques dudit objet, comportant une 

5 selection dans une premiere image st£r6otaxique (10) d'une premiere 

region d'interet (5d) dite "cible", ainsi qu'un appariement de ladite 
premiere region avec une deuxieme region d'interet (5g) homologue de 
la premiere et apparaissant dans une deuxieme image st6r6otaxique 
(12), caracteris6 par le fait que, les images stereotaxiques etant 

10 numerisees, on selectionne (400) dans 1'etape de selection, un pixel 

cible (60d) dans ladite region d'interet cible (5d), et dans 1'etape 
d'appariement, on genere (401), autour du pixel cible selectionne, une 
fenetre cible (6 Id) de caracteristiques dimensionnelles choisies et 
contenant ladite region d'interet cible, on determine (402) dans la 

15 deuxieme image, selon un critere de selection predetermine, un 

ensemble de pixels, on genere (403) autour de chaque pixel 
selectionne une deuxieme fenStre de memes caracteristiques 
dimensionnelles que ladite fenetre cible, on effectue un traitement de 
correlation (404) entre les niveaux de gris des pixels de chaque 

20 deuxieme fenetre et les niveaux de gris des pixels de la fenetre cible, 

de fa?on a obtenir une valeur de correlation pour chaque deuxieme 
fenetre, et on identifie la region d'interet homologue de la region 
d'interet cible a partir de l'analyse de I'ensemble des valeurs de 
correlation ainsi obtenues, de fa$on a minimiser les risques d'erreur 

25 d'appariement entre les regions d'interet homologues. 

2. Procdde selon la revendication 1, caractdrise par le fait que 
Tanalyse des valeurs de correlation obtenues comporte la selection 
d un nombre choisi de maxima ou de minima de correlation, la 
selection de la region d'interet homologue s'effectuant parmi celles 

30 dont la valeur de correlation associee est 1'un de ces maxima ou 

minima de correlation. 

3. Procede selon la revendication 2, caracterise par le fait que 
l'analyse des valeurs de correlation obtenues comporte en outre une 
determination de la dynamique des maxima ou minima seiectionnes, et 
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la comparaison des valeurs de dynamique obtenues a un seuil. 

4. Proc£de selon Tune des revendications 1 a 3, caracterise par le 
fait que le traitement de correlation comporte un traitement de 
correlation normee. 

5. Proc6de selon Tune des revendications 1 h 3, caracterise par le 
fait que le traitement de correlation comporte un traitement de 
difference normee. 

6. Proc6de de localisation d'un element d'interet contenu dans un 
objet tridimensionnel a partir des positions de regions d'interets 
homologues correspondant audit Element d'interet et apparaissant dans 
un ensemble damages st6reotaxiques dudit objet, comportant une 
selection dans une premiere image stereotaxique (10) d'une premiere 
region d'interet (5d) dite "cible", ainsi qu'un appariement de ladite 
premiere region avec une deuxieme region d'interet (5g) homologue de 
la premiere et apparaissant dans une deuxieme image stereotaxique 
(12), caracterise par le fait que, les images stereotaxiques etant 
numerisees, on s61ectionne(400) dans retape de selection, un pixel 
cible dans ladite region d'interet cible, et dans retape d'appariement, 
on genere (401), autour du pixel cible selectionne, une fenStre cible 
(6 Id) de caracteristiques dimensionnelles choisies et contenant ladite 
region d'interet cible, on determine (409), pour ladite region d'interet 
cible, a partir des valeurs de niveaux de gris des pixels de ladite 
fenetre cible, un premier ensemble de valeurs numdriques 
respectivement representatives de caracteristiques predetenmindes 
caracterisant toute region d'interet, et on stocke ce premier ensemble 
de valeurs numeriques, on determine (402) dans la deuxieme image, 
selon un critere de selection predetermine, un ensemble de pixels 
appartenant respectivement h des regions d'interet dites "candidates" 
apparaissant dans la deuxieme image, et on genere (403) autour de 
chacun de ces pixels selectionnds une deuxieme fenetre contenant 
ladite region d'interet candidate correspondante, on determine (410), 
pour chaque region d'interet candidate, a partir des valeurs de niveaux 
de gris des pixels de la deuxieme fenetre associee, un deuxieme 
ensemble de valeurs numeriques respectivement representatives 
desdites caracteristiques predeterminees, on stocke ce deuxifeme 
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ensemble de valeurs numeriques, et on identifie la region d'int6ret 
homologue de la region d'interet cible a partir d'un traitement de 
comparaison (411) entre ledit premier ensemble de valeurs numeriques 
et chaque deuxieme ensemble de valeurs numeriques, de fa9on a 
minimiser les risques d'erreur d'appariement entre les regions d'interet 
homologues. 

7. Procede selon la revendication 6, caracterise par le fait que 
lesdites caracteristiques predeterminees comprennent des 
caracteristiques de forme, de contraste, de gradient. 

8. Procede selon la revendication 6 ou 7, caracterise par le fait 
que le traitement de comparaison comporte un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

9. Procede selon Tune des revendications precedentes, caracterise 
par le fait qu'on determine dans la deuxieme image stereotaxique une 
zone epipolaire contenant au moins le segment epipolaire relatif au 
pixel cible, et par le fait que les pixels de la zone epipolaire 
contiennent les pixels selectionnes. 

10. Procede selon Tune des revendications precedentes, 
caracterise par le fait que lesdits pixels selectionnes sont choisis 
parmi les pixels de la deuxieme image presentant des maxima ou 
minima de niveaux de gris dont la dynamique est superieure a un seuil 
predetermine. 

11. Procede selon Tune des revendications precedentes, 
caracterise par le fait que l'appariement comporte prdalablement un 
filtrage des images stereotaxiques, en particulier du type dit "chapeau 
haut de forme". 

12. Procede de localisation d'un element d'interSt contenu dans un 
objet tridimensionnel h partir des positions de regions d f interets 
homologues correspondant audit element d'interet et apparaissant dans 
un ensemble damages stereotaxiques dudit objet, caracterise par le fait 
que, les images stereotaxiques etant numerisees, il comprend une 
selection dans une premiere image st6reotaxique (10) d'une premiere 
region d'interet (5d) dite "cible", une selection (405), dans une 
deuxieme image stereotaxique (12), a partir d f un premier appariement 
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automatique, d'au moins une deuxifeme region d'interet (5g) suceptible 
d'etre homologue de la region cible, la d6termination de la position 
spatiale d'un pixel candidat d'un element d'interet candidat 
correspondant h ces deux regions d'interet, la determination (406) de 
la position dans une troisieme image st<5r£otaxique, d'un pixel projet€ 
correspondant a la projection dans cette troisieme image sterdotaxique 
du pixel candidat, et un deuxieme appariement automatique entre la 
region d'interet cible et un voisinage du pixel projet6 d£finissant une 
region d'interet projetee, de fa9on a minimiser les risques d'erreur 
d'appariement entre les regions d'interet homologues. 

13. Proced6 selon la revendication 12, caract£ris6 par le fait qu'on 
selectionne dans la deuxi&me image st6reotaxique, k partir dudit 
premier appariement automatique, un ensemble de deuxiemes regions 
susceptibles d'etre homologues de la region cible, on dftermine la 
position spatiale de chaque pixel candidat correspondant & chaque 
paire de regions d'interet form6e par la region cible et Tune des 
deuxiemes regions, ainsi que la position dans la troisieme image de 
chaque pixel projete correspondant, et on effectue le deuxieme 
appariement automatique entre la region cible et un voisinage de 
chaque pixel projet6. 

14. Procede selon la revendication 13, caracterise par le fait qu'on 
ordonne l'ensemble des deuxiemes regions obtenues dans l'ordre 
decroissant de leur probability a etre la region homologue de la region 
cible, on ordonne l'ensemble des regions projetees obtenues dans 
l'ordre decroissant de leur probability & 6tre la region homologue de la 
region cible, et on sdlectionne la region homologue de la region cible 
comme 6tant celle dont le produit des rangs respectifs dans les deux 
ordonnancements est minimum. 

15. Proc6d6 selon Tune des revendications 12 & 14, caract£ris6 par 
le fait que le premier appariement automatique est un appariement tel 
que defini dans Tune des revendications 1 & 11. 

16. Troc6d6 selon Tune des revendications 12 k 15, caxacteris6 par 
le fait que le deuxifeme appariement automatique comporte un 
traitement de correlation tel que d6fini dans Tune des revendications 1 
h 5, entre le voisinage de chaque pixel projetd et une fenetre cible 
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contenant la r6gion d'interet cible, ou un traitement de comparaison, 
tel que d£fini dans Tune des revendications 6 a 11 entre la region 
d'interet cible et chaque region d'interet projetee. 
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